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Current-voltage (I-V) characteristics of quench condensed, 
superconducting, ultrathin Bi films in a magnetic field are re- 
ported. These I-V's show hysteresis for all films, grown both 
with and without thin Ge underlayers. Films on Ge underlay- 
ers, close to superconductor- insulator transition (SIT), show 
a peak in the critical current, indicating a structural transfor- 
mation of the vortex solid (VS). These underlayers, used to 
make the films more homogeneous, are found to be more ef- 
fective in pinning the vortices. The upper critical fields (Bc2) 
of these films are determined from the resistive transitions in 
perpendicular magnetic field. The temperature dependence 
of the upper critical field is found to differ significantly from 
Ginzburg-Landau theory, after modifications for disorder. 
PACS Numbers : 73.50.-h, 74.40. +k, 74.80.Bj 



I. INTRODUCTION 

Transport properties of disordered ultrathin films have 
been studied extensively over the last decade especially in 
the context of disorder driven and magnetic field driven 
superconductor-insulator transition (SIT) This tran- 
sition is found to occur at a particular value of disorder 
corresponding to a critical resistance Rc, which clusters 
close to ^ h/4e^(the quantum of resistance for Cooper 
pairs), for both field driven and disorder driven transi- 
tions. Earlier studies on Bi, quench-condensed on under- 
layers of Ge, showed a "homogenous" type of SIT which 
was characterized by a strong suppression of the critical 
temperature, Tc as disorder was increased. The presence 
of a thin (< 10 A) underlayer of Ge, is conventionally 
thought to improve the wetting properties of the film, 
and thereby assist homogeneous growth Q|. However, 
several studies have indicated that the underlayer is not 
inert, as conventionally assumed, and may actually play 
an active role in determining the transport properties of 
the film Though most studies have concentrated on 
the transition itself, little effort has been made to under- 
stand the physics of superconductivity, such as vortex 
dynamics, in ultrathin films with different underlayers, 
which differs significantly from its manifestation in bulk 
material. In this work, we demonstrate some of the dif- 
ferences. In a previous publication, we have investia- 
gated several underlayers for these films including solid 
Xenon to investigate screening effects. 

In this paper we present observations of some anoma- 
lies in the I-V's in a magnetic field and the Bc2{T) of 
Bi films, which are in the vicinity of the disorder driven 
SIT. This work compares films without and with a Ge 



underlayer. From current-voltage (I-V) characteristics in 
the superconducting state, measured in finite magnetic 
field, it is possible to infer the dynamics of vortices. 



II. EXPERIMENTAL RESULTS 

Experiments were carried out on films near the dis- 
order driven SIT, with and without Ge underlayers si- 
multaneously to determine the effect of underlayers on 
the vortex motion in these films. Upper critical fields 
Bc2{T) of these films were also measured, from the re- 
sistive transitions in perpendicular magnetic field, and 
their temperature dependance was found to differ signif- 
icantly from Ginzburg-Landau theory Q. The experi- 
ments were done in a UHV cryostat, custom designed 
for in-situ experiments, described in The cryostat is 
ptimped by a turbomolecular pump backed by an oil-free 
diaphragm pump. A completely hydrocarbon free vac- 
uum < 10^^ Torr can be attained. The substrates are 
amorphous quartz of size 2.5cm X 2.5cm and mounted 
on a copper cold finger whose temperature can be var- 
ied down to 1.8K by pumping on the liquid helium bath. 
The material {Bi) is evaporated from a Knudsen-cell with 
a pyrolytic Boron Nitride crucible, of the type used in 
Molecular Beam Epitaxy (MBE). Bi is evaporated from 
the cell at 650° C, into a 4-probe resistivity measurement 
pattern defined by a metal mask in front of the substrate. 
Successive liquid helium and liquid nitrogen cooled jack- 
ets surrounding the substrate reduce the heat load on 
the substrate, and provide cryo-pumping. This produces 
an ultimate pressure ~ 10^^" Torr in the system. The 
fiux reaching the substrate is controlled using a carefully 
aligned mechanical shutter in the nitrogen shield. The 
thickness of the film is increased by small amounts, by 
opening the shutter for a time interval corresponding to 
the desired increase in thickness. A quartz crystal thick- 
ness monitor measures the nominal thickness of the film. 
Electrical contacts to the film are provided through pre- 
deposited platinum contact pads (^ 50 A thick). Ge un- 
derlayers are deposited over a region of the substrate (a- 
quartz) before loading the substrate into the cryostat. 
Separate electrical connections to films on Ge and on 
bare a-quartz allow us to study both the films simulta- 
neously. I-V's and electrical resistance measurements are 
done using a standard d.c. current source (Keithley model 
220/224) and nanovoltmetcr (Keithley model 182) and 
elctrometer (Keithley model 6514). 

In this work, we have measured the critical current Ic 
for different films at T = 2. OK. We have evaluated the 
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power dissipated by the current in the superconducting 
films to be a few hundreds of nanowatts at the highest 
currents used (i.e. near the critical current). The cryo- 
stat has been designed to have a cooling power of at least 
1 mW at this temperature. The substrates on which the 
films are deposited, are mounted on a copper coldfingcr 
with silver loaded Apiezon grease, ensuring a low ther- 
mal resistance. Under such circumstances, we believe the 
differences in the cooling power and ohmic heating to be 
large enough to ensure that the temperature fluctuations 
due to ohmic dissipation can be safely ignored. 

When the normal state sheet resistance, of the film 
is close to h/4e^, the combined effect of quantum fluctu- 
ations and disorder prevents the formation of a vortex 
solid (VS) phase, and such an absence of a VS phase 
has recently been reported in Bi films close to SIT 
However, we find that this type of behavior can be modi- 
fied when the underlayer is changed. The purpose of this 
paper is to present our results which establish the influ- 
ence of the underlayer on the vortex phase, and others, 
which go beyond those reported in ref. [6]. Specifically, 
these are the observation of a structural phase transition 
in the vortex phase, which manifests itself as a 'peak ef- 
fect', and the significant deviations in the critical fields 
of these ultrathin superconducting films, when compared 
to Ginzburg-Landau theory. Keeping this in mind, we 
study I-V characteristics in magnetic field to understand 
the effect of magnetic field on the vortices and how it 
differs for films on different substrates. In the following, 
after suitable preliminaries, we present our results and 
discuss their implications. 

The I-V curves are hysteretic as shown in Fig. 1. i.e. 
when the current is increased from the zero value, the 
resistance jumps to the normal state value at a critical 
current, Ic. Upon reducing the current from the nor- 
mal state, the voltage returns to zero not at Ic, but at 
a lower value. This behavior is reminescent of a under- 
damped Josephson junction, where the intrinsic junction 
capacitance, and the finite leakage resistance to quasi- 
particles can be modeled as shunts to an ideal junction. 
The lower current in this case is the retrapping current. 
This suggests that our films can be considered as an array 
of Josephson junctions, which are shunted by an equiv- 
alent resistance and the equivalent intrinsic junction ca- 
pacitance. Consequently the resistively and capacitively 
shunted junction (RCSJ) model has be used to de- 
scribe the hysteretic behavior |^ in an earlier publica- 
tion. 

Fig. 2 shows the I-V characteristics for a 26 A Bi film 
on (a) on bare a-quartz substrates and (b) 10 A Ge un- 
derlayer, both on a linear scale, at low voltages. In this 
figure only the positive current sweep is shown for clarity, 
though all the curves show hysteretic behavior. All the 
I-Vs are recorded at T = 2. OK. The differences are de- 
scribed below. In (a), at zero field, one can clearly see the 
dissipationless supercurrent region. Upon increasing the 



field, the system first becomes dissipative due to depin- 
ning of the vortices, caused by the Lorentz force acting 
on them due to the current Idep, although not yet nor- 
mal. At this point the TV becomes nonlinear, and the 
departure from linearity is clearly seen in this figure. We 
note here that the definition of the depinning current is 
consistent with the one conventionally used in the litera- 
ture . For a film of identical thickness, the Ic is higher 
on the Be underlayer, as can be seen from Fig. 2(b). We 
also note that the I-V is non-linear in the vicinity of the 
Ic. The non-linearity attributable to depinning seems to 
set in at lower currents for the films on Ge. However, 
the Ic for films on Ge is higher. In both these films, the 
dissipation is finite at any applied field, as evidenced by 
a linear regime of the I-V, from which a resistance can 
be discerned. This is possibly due to the relatively high 
temperature of 2 K, when compared with published liter- 
ature . This we attribute to thermally activated vortex 
motion. 

Fig. 3 (a) and (b) show the I-V's on a semi-log scale, 
until the normal state resistance is reached, where all the 
curves at different magnetic fields meet. For the films on 
bare substrates, only one critical current can be unam- 
biguously identified. For the films on Ge underlayers, two 
different step like increases in the voltage can be clearly 
identifed for three different field values. It is only after 
the second step that the film resistance becomes the nor- 
mal state resistance, and therefore this should be defined 
as the critical current. Hence the first step has to be 
associated with some sort of structural transition in the 
vortex solid. Recall that the depinning current is lower 
than this current at which the first step in the voltage 
appears. 

Fig. 4 shows the critical current as a function of nor- 
malized temperature, at zero field for films on Ge as 
well as on bare quartz. A peak in the critical current in 
the vicinity of the transition temperature is clearly vis- 
ible, for the films on Ge. An identical phenomenon has 
been extensively invesitaged in layered superconductors 
such as 2H-NbSe2 |^ and several other materials such as 
amorphous NbsGe and MosSi films, jl^] This has been 
given the name "peak effect". The pinning force density 
Fp exhibits a prominent peak slightly below the upper 
critical field, Bc2 or near the critical temperature. This 
effect is attributed to elastic instabilities generated by 
local fluctuations of the pinning forces, which induce a 
rapid softening of the flux line solid. This softening al- 
lows the flux lines to conform readily to a configuration 
that locks it to the inhomogeneities, so that the critical 
current Ic increases [|ll|. So the observation of a peak 
in Ic implies that an elastic flux line solid transforms to 
a plastic solid in the vicinity of the peak. For the films 
on Ge, we identify the first transition with the "peak ef- 
fect" . On the other hand, for films on bare substrate at B 
= G, the system shows the usual depinning transition 
with some finite dissipation and non-linear I-V's, until Ic, 
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where it becomes normal. In our films, the critical fields 
are higher than what we can attain using our magnet. 
Hence we have been unable to observe the peak effect in 
the I-H plane. However, the observation in the I-T plane 
is unambiguous proof of the peak effect |Q . 

There is another important difference between the 
films on Ge and bare substrates. The Ic when plotted 
against the applied field B, shows a maximum for films 
on bare substrates as shown in Fig. 5. This effect is 
therefore real, as the peak is much larger than the er- 
ror bars. Such an effect has been observed in Josephson 
junction arrays, where the cirtical current shows a depen- 
dence on the field for small magnetic fields . In such J J 
array structures, a relevant variable is the number of free 
vortices per unit cell, and is denoted as f (/ = Bao^/^o), 
where B is the applied field. Since we have sufficient 
reason to believe that our films can be understood as a 
random J J array as discussed above this rea- 

sonable explanation for this effect. From the measured 
Ic, and using the Ambegaokar-Baratoff relation |l5), we 
can estimate the number of junctions in our film to be 
a million. Knowledge of the applied field, and using our 
published RHEED data (5[|l|,0 to estimate a^ (of the 
order of 2000 A), gives us a value for f of 0.4, which is 
consistent with the values reported in the literature. It 
is intriguing that such a frustration effect is not observed 
for films on Ge underlayers. We do not know the reason 
for this at this time. We are however in the process of 
investigating other films and underlayers in order to aid 
our understanding of this phenomenon. 

We now discuss the critical fields of these films. The 
temperature dependence of the upper critical field, Bc2 
of a typical type II superconductor can be shown, by 
Ginzburg-Landau theory to be of the form pO[|, 



Se2 =Sc2(0)(l-t2) 



(1) 



where Bc2(0) is the upper critical field at T=0 and t (= 
T/Tc) is the reduced temperature. Close to the tran- 
sition temperature (t ~ 1), the relation above can be 
approximated by (1 — i). 

Fig. 6 shows a representative Bc2 vs t plot for a 26 A 
film where the data points are indicated by triangles, and 
the solid line is a fit discussed below. The data show a 
curvature opposite to that expected by GL theory near 
Tc- The data is found to fit the following equation. 



Bc2 = B,2m^-tT 



(2) 



where a = 1.175. Even at t ~ 1, we do not observe the 
predicted linear dependence. In Fig. 6, the solid line is 
a fit to equation 3 with a = 1.175, the long-dashed line 
is the GL expression, medium-dashed line is the linear 
dependence, and the short-dashed line is a = 1.45. This 
last value of a is found to characterize the melting line 
in extreme type II superconductors, such as the high Tc 
materials [p3. All the films studied (thickness ranging 



from 26 to 105 A) show the same curvature with a 
1.14 ± 0.05. 



III. DISCUSSION 

The transport properties in films such as the ones stud- 
ied here are examples of transport in heterogeneous me- 
dia |16-18|, i.e. as the problems of dynamics of percola- 



tion networks. The current path in such films is clearly 
percolative, and one cannot rule out such effects in the 
motion of vortices. A possible model for these films seems 
to be a random composite network (RCN) , in which both 
superconducting and normal regions exist |l9| ] . Attempts 
have been made to quantitatively interpret the I-V's in 
such structures [|l8|. Currently we are in the process of 
incorporating the two dimensional nature of our system, 
and the effects of a magnetic field to quantitatively un- 
derstand our I-V's. 

The I-V's are sensitive to the deposition temperature. 
All the TV's discussed above were obtained from films for 
which the temperature of deposition was 15 K. Several 
dramatic differences are seen in the TV's of films that 
are deposited at 4 K or 1.5 K. These results along with 
a quantitative model for the I-V's will be addressed in a 
separate publication. 

For a disordered superconductor, the temperature de- 
pendence of the critical field gets modified due to the 
interplay between pinning, and thermal fluctuations. In 
a weakly disordered superconductor, the temperature de- 
pendence of the upper critical field shows significant de- 
viations from this relation only at low temperatures. The 
interplay between pinning and fluctuations increases the 
critical flelds. Such increases have been successfully ex- 
plained by the Werthamer-Helfand-Hohenberg (WHH) 
theory ||2l| , which is based on the Ginzburg-Landau for- 
malism. However, the WHH model has not been success- 
ful in explaining the Bc2(r) of amorphous transition 
metal alloys which often show an upward curvature. 

Later models taking into consideration, the effect of 
localization and electron-electron interactions ||2^, were 
successful in explaining the upward curvature of the Bc2 
at lower temperatures. We have determined Bc2(T) of 
our fllms, from the resistive transitions in a perpendicu- 
lar magnetic field. The convention that we have followed 
is to define Bc2 as the field at which the sheet resistance 
is half its normal state value, R^, which is the generally 
accepted convention. The current values used to record 
the resistive transition of the films are small few nano 
amperes) and the magnetic field produced by them is ex- 
tremely small (milligauss) compared to the critical fields 
(Tesla) of these films. 

Although the theories discussed above [ pT|j23| explain 
the upward curvature of Bc2 at lower temperatures very 
well, they predict only a linear dependence (i.e. (1 — t)) 
of Bc2 near Tc. None of them has any explanation for 
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this peculiar behavior observed in our films. A power of 
^ 1.14, which is between the value of 1.45 for a melting 
line and 1.0 for GL theory suggests that a new expla- 
nation may be required to understand this data. The 
Pauli limiting field or the Clogston-Chandrasekhar limit 
for disordered superconductors is given by Hp = 1.84 Tc 
Tesla The value of Hp is approximately 5.3 T even 
the for the thinnest films studied here. In our experi- 
ments we have not crossed that field at any time. 

IV. CONCLUSIONS 

In conclusion, the I-V's of quench condensed Bi films 
with and without underlayers of Ge show evidence for 
very interesting vortex dynamics. In particular, we have 
presented data that show a peak effect, which has been 
observed for the first time in such ultrathin films. We 
also show effects that can be attributed to frustration in 
a random JJ array. These effects are sensitive to the mi- 
crostructure of the film, and the underlayers used to de- 
posit the film. The upper critical fields are found to have 
a power law temperature dependence with a power of ~ 
1.14, which needs further investigation. This, together 
with the peak effect and frustration effects, suggests a 
novel behavior in such ultrathin disordered films, which 
has not been addressed. More experiments on highly dis- 
ordered films with Rj^ h/4e^ are needed to completely 
understand the vortex dynamics. 
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FIGURE CAPTIONS 

(1) Fig. 1. I-V characteristics of a 26 A film on Ge 
underlayer showing hysteresis for (1) zero field (2) 364 G 
and (3) 7981G. The curves were recorded at T = 2.0K. 

(2) Fig. 2. I-V characteristics of 26 A Bi film on (a) 

bare quartz and (b) on a 10 A Ge underlayer on a lin- 
ear scale at low voltages. All I-Vs were recorded at T = 
2. OK. Different curves indicate different field values. 

(3) Fig. 3. I-V characteristics of 26 A Bi film on (a) 
bare quartz and (b) on a 10 A Ge underlayer on a semi- 
log scale, taken until the films are driven normal by te 
current. All 1-Vs were recorded at T = 2. OK. Different 
curves indicate different field values. 

(4) Fig. 4. Temperature dependance of the critical 
current (Ic) for films on bare substrates and films on Ge 
under layers. 

(5) Fig. 5. Field dependance of the critical current 
(Ic) for the film of Fig. 2 (b). 

(6) Fig. 6. Bc2 vs t for a 26 A film. The solid line is the 
fit to equation 3 with a = 1.175. See text for discussion. 
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